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1  INTRODUCTION 


1  Introduction 


Modeling  of  high-enthalpy  plasmas  is  essential  to  the  fundamental  understanding  and  de¬ 
velopment  of  many  applications  of  interest  to  the  Air  Force,  including  aerodynamic  flow 
control,  re-entry  and  hypersonic  air-breathing  vehicles,  as  well  as  plasma  discharges  for 
electric  propulsion  devices.  As  applied  to  such  cases,  modeling  can  aid  in  improved  aerody¬ 
namic  performance,  thermal  protection  systems,  combustion  enhancement,  and  thruster 
efficiency  among  others.  The  difficulty  in  modeling  such  plasmas  is  the  high-degree  of 
nonequilibrium  between  the  various  plasma  components.  Energy  exchange  rates  between 
the  various  internal  modes,  for  example,  may  be  comparable  to  the  transport  time  scales, 
necessitating  multi-temperature  and  collisional-radiative  (CR)  models.  Furthermore,  mo¬ 
mentum  exchange  between  the  various  components  may  be  slow  requiring  a  multi-fluid 
description  of  the  plasma.  These  complexities  are  further  compounded  when  a  continuum 
description  of  the  plasma  is  no  longer  valid,  making  it  necessary  to  augment  traditional 
CFD  approaches  with  kinetic  or  moment  equations. 


In  a  given  flow  field  simulation,  the  plasma  conditions  can  vary  greatly,  making  it  neces¬ 
sary  to  combine  multiple  descriptions  of  the  plasma  as  determined  by  the  local  conditions, 
and  in  particular,  the  degree  of  nonequilibrium.  For  example,  the  plasma  transport  may 
transition  between  rarefied  and  continuum  flow,  requiring  appropriate  models  for  each  case 
through  hybrid  techniques.  Such  hybridization  can  also  be  important  to  the  chemistry  of 
the  flow.  In  the  case  of  electronic  excitation  and  ionization,  detailed  collisional-radiative 
models  may  be  essential  in  accurately  capturing  the  flow  in  certain  regions,  while  un¬ 
necessary  and  otherwise  computationally-burdensome  in  regions  in  which  nonequilibrium 
effects  become  negligible.  In  such  cases,  a  hybridization  approach  between  a  CR  model 
along  with  a  simple  redistribution  of  levels  according  to  a  Boltzmann  distribution  where 
appropriate  can  yield  significant  performance  gains. 


In  particular,  it  is  desired  to  develop  models  based  on  a  multi-temperature,  multi-component 
description  of  the  plasma  combined  with  collisional-radiative  models  in  order  to  describe 
the  chemical  kinetics.  Recent  research  on  shocktube  simulations  in  monatomic  gases  has 
provided  valuable  insight  into  the  kinetics  and  dynamics  of  the  collisionally-dominated 
regime  of  the  plasma.  However,  it  is  also  important  to  investigate  regimes  in  which  ra¬ 
diative  process  can  dominate  collisional  processes.  One  such  example  is  where  radiative 
recombination  begins  to  dominate  three-body  recombination  when  there  is  a  lack  of  third 
body  collision  partners. 


The  current  work  is  focused  on  the  development  of  CR  models  for  monatomic  gases.  The 
models  must  be  detailed  enough  to  resolve  the  degree  of  nonequilibrium  in  the  plasma 
while  remaining  robust  and  efficient  enough  to  permit  coupling  with  unsteady,  multi¬ 
dimensional  flow  solvers.  We  first  develop  the  models  for  argon,  krypton  and  xenon  and 
then  validate  them  using  shocktube  data. 
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2  General  model 

We  begin  by  describing  a  general  CR  model  which  is  then  applied  specifically  to  Ar,  Kr 
and  Xe.  Since  detailed  cross  sections  for  all  transitions  are  not  known  by  either  theoretical 
or  experimental  means,  generalized  models  can  provide  useful  approximations  to  otherwise 
unknown  quantities.  These  semi-empirical  models  are  validated  on  transitions  where  data 
are  available  and  can  thus  provide  trends  and  insights. 

The  processes  considered  can  be  represented  by  the  following  equations, 


Ar(i)  +  e  Ar (J)  +  e  (1) 

Fji 

Ar(i)  +  Ar(l)  Ar  (j)  +  Ar(l)  (2) 

Lji 

Ar (i)  +  e~  S=^  Ar+  +  e_  +  e_  (3) 

Oi 

Ar(i)  +  Ar(l)  ;=A  Ar+  +  e_  +  Ar(l)  (4) 

Wi 

Ar (i)  +  hu^  v  x  Ar(j)  (5) 

Aj  i  Aj  i 

Ar(z)  +  hvi  s  x  Ar+  +  e_  (6) 

A  %Ri 


where  the  rate  coefficients  defined  by  [33]  are  summarized  in  Table  1.  Although  the 
processes  in  this  section  are  expressed  for  argon,  they  are  equally  applicable  for  krypton 
and  xenon.  The  convention  adopted  here  is  j  >  i  such  that  Cl3  and  KtJ  represent  the 
rates  for  excitation  from  level  i  to  j  while  F]t,  Lji,  and  A:]t  are  the  rates  for  de-excitation 
from  level  j  to  i.  The  levels  considered  for  the  models  presented  herein  include  only  those 
contained  in  the  lowest  lying  manifolds,  beyond  which  the  ionization  limit  is  approximately 
1  eV  or  less.  Thus  the  lower- lying  states  considered  here  make-up  the  bottleneck  to  the 
ionization  regime  and  contain  most  of  the  energy  of  bound  states.  Note  also  that  neutral 
atoms  have  two  ionization  potentials  owing  to  the  fact  that  the  levels  of  argon  are  split 
according  to  two  possible  configurations  for  the  core  angular  momentum  jc.  As  ionization 
and  recombination  can  proceed  only  between  levels  having  the  same  value  of  jc,  the  two 
lowest  lying  levels  of  the  corresponding  ions  must  be  accounted  for  individually.  Thus  in 
the  current  model,  these  two  levels  are  treated  as  separate  species  and  convected  as  such. 

2.1  Atom  impact  processes 

Heavy-particle  impact  cross  sections  are  scarce  in  the  literature,  requiring  us  to  rely  upon 
semi-empirical  models  including  that  of  Drawin  [19] .  However,  caution  must  be  exercised 
when  using  such  models  as  are  known  to  over-estimate  the  impact  cross  sections  by  at  least 
an  order  of  magnitude  [10].  In  several  applications  of  general  interest,  including  shock- 
heated  plasma,  atom-atom  collisions  are  the  only  means  for  initial  electron  production  and 
yield  the  priming  electrons  that  eventually  trigger  inelastic  electron  collisions.  Excitation 
from  ground  state  in  particular  plays  a  significant  role  in  determining  the  time  necessary 
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Rate 

Coefficient 

Process 

Cn 

collisional  excitation  by  electrons 

Kij 

collisional  excitation  by  ground  state  atoms 

Fji 

collisional  de-excitation  by  electrons 

Lji 

collisional  de-excitation  by  ground  state  atoms 

Si, Vi 

collisional  ionization 

Oi,Wi 

three-body  recombination 

Ri 

radiative  recombination 

Aij 

transition  probability/spontaneous  emission  (Einstein  coefficient) 

Fij 

bound-bound  optical  escape  factor 

K 

bound-free  optical  escape  factor 

k  ■ 

^ei 

electron-ion  collisions 

h 

r^en 

electron-neutral  collisions 

^brems 

bremsstrahlung  emission 

Table  1:  Summary  of  rate  coefficients  for  collisional-radiative  model. 


for  the  onset  of  an  electron  avalanche  behind  a  viscous  shock  for  example.  Starting  with 
Drawin’s  formula 


<7'(£)  =  4?rao 


'if 


mAr  d2  f  2 me 

S  Jij  , 
mu  m-Ar  +  "A 


e/Ae. 


1  + 


2m  e 


mAr+me 


(e/Ae. 


(7) 


where  /#  is  the  ionization  potential  of  the  hydrogen  atom,  aQ  is  the  Bohr  radius  and  £  =  6 
is  the  number  of  optical  electrons  of  argon.  For  the  energy  ranges  under  consideration, 
Eq.  (7)  is  well-approximated  by  a  linear  function, 


where 


°ij(£)  =  Pij(£  ~  A%), 

(8) 

a*  a„„2{!h)2  C2f  2me 

Pi j  —  Fxa0  3  s  Jij 

J  efj  rrxH 

(9) 

and  fij  is  the  oscillator  strength  of  the  transition. 

Less  sensitive  are  the  ionization  cross  sections.  The  formula  of  Drawin  (Eq.  (7),  replacing 
fij  =  1  and  £ij  =  I,,  the  ionization  potential)  can  be  used  to  obtain  the  cross  sections  for 
ionization  for  the  excited  levels1, 


of(e)  =  4t ra„ 


2 


TTlAr  ^ 2  2me 

rriff  mAr  +  ure 


ejU  -  1 


2m  e 


mAr+me 


(V-fi  -  d) 


2  ‘ 


1for  further  details,  the  reader  is  directed  to  [4]  and  [33]  and  the  references  therein 


(10) 


Final  Report 


-  5  - 


FA  8655-10-1-3076 


2  GENERAL  MODEL 


2.2  Electron  impact  excitation 


Zatsarinny  &  Bartschat  [35]  have  published  cross  sections  using  R-matrix  calculations 
for  most  of  the  noble  gases  which  agree  favorably  with  the  detailed  experimental  results 
of  Buckmann  et  al  [3].  Recent  experiments  by  have  provided  further  validation  of  the 
R-matrix  calculations.  For  transitions  for  which  such  detailed  computations  are  not  yet 
available,  the  formula  of  Drawin, 


otj(e)=4*a^ 


-1) 


u, 


( afjfij  [lHj £ij)2  ln(1.25 Ajt%) 

\< 

[afjiuij  +  1  )/u% 


may  be  used  with  atf-fij  —  1,  aT  —  1,  and  cxfj  —  1. 


(11) 


2.3  Electron  impact  ionization 


For  electron  impact  cross  sections,  we  rely  upon  the  Deutsch-Mark  formalism  [6]  which 
combines  a  semi-empirical  Gryzinski-type  energy  dependence  of  the  cross  section  with 
quantum  mechanically  calculated  molecular  structure  information.  The  cross  sections 
can  be  expressed  as 

<7i{e)  =  gntnrleZrufie),  (12) 

where  gne  are  the  reduced  weighting  factors,  rn£  are  the  radii  of  the  valence  electron  (n,  £ 
are  the  principal  and  orbital  quantum  numbers  describing  the  electronic  state),  and 


fi(e)  =  d— 
£ 


(ejh  ~  1 
\e/Ii  +  1 


a 


b  +  c 


|)ln(2.7+(£//,-l)1/2) 


(13) 


the 


2.4  Elastic  collisions 

Mitroy  [25]  used  modified  effective  range  theory  (MERT)  to  develop  a  semi-empirical 
fit  for  momentum  transfer  cross  sections  involving  noble  gases,  both  pure  as  well  as  in 


Atom 

Dipole,  ad 

Polarizability 

[Oq]  Quadrupole,  aq  [a®] 

Ne 

2.69 

7.52 

Ar 

11.07 

52.25 

Kr 

17.06 

97.39 

Xe 

27.66 

209.85 

Table  2:  Dipole  and  quadrupole  polarizabilities  of  the  noble  gases  ([28]). 
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mixtures  with  H2.  According  to  the  theory,  the  cross  section  can  be  expressed  in  terms 


of  the  phase  shifts  for  the  partial  waves  which  are  given  by 

tan  770  =  —Ak[  1  +  (4^/3) k2  In  k]  —  ( iraci/3)k 2  +  Dk 3  +  Fk4  (14) 

tan  771  =  aioidk2  —  Aik 3  +  {b\o?d  +  c\a2)k4  +  Hk 5  (15) 

tan 772  =  a2adk 2  -  ( b2a2d  +  c2a2)k4  +  A2k5  (16) 

ta.nr]e  =  aiadk2  +  (bea2  + C£aq)k4  (17) 

(18) 


where  £  is  the  electron  angular  momentum  number.  The  coefficients  ae,  bp  and  cp  are  given 
by 


7T 

(2f  +  3)(2£  +  1)(2£  —  1) 

tt[15(2£  +  l)4  -  140(21?  +  l)2  +  128] 
[{2i  +  3)(2£  +  1)(2£  -  1)]3(2£  +  5)(2£  -  3) 
3  d£ 

(2£  +  5)(2£-3)' 


(19) 

(20) 

(21) 

(22) 


In  terms  of  the  phase  shifts,  the  momentum  transfer  cross  section  is  given  by 

<Tm  =  'P'X^+1  )sin2(%  -  %+i)-  (23) 

The  dipole  ad  and  quadrupole  aq  polarizabilities  for  the  noble  gases  are  summarized  in 
Table  2.4. 


2.5  Photorecombination/Photoionization 


We  start  with  photorecombination,  which  can  play  an  important  role  in  radiative  cool¬ 
ing.  The  cross  section  for  photorecombination  can  be  found  from  the  cross  section  for 
photoionization  under  equilibrium  conditions  (Zel’dovich  &  Raizer  [37]), 


<7c,i(U)  = 


h2v2 


g+  mlv^c2 


(24) 


Utilizing  the  relation  hu  =  mev2  / 2  +  £,  =  £  +  £,,  the  cross  section  associated  with  the 
ground  state  is  given  by 


<7„,  i(e) 


9i  (s  +  £i)2 
g+  2  emec2 


x 


{3.5  x  10"21 
2.8  x  lO"20 


lH 


£  ~\~  £  1 


0  <  c  <  2IH  —  £\ 


£  ^ >  2IH  —  £\ 


(25) 
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while  all  others  are  computed  using 


=  7* 


9i  (g  +  O)2 
g+  2  emec2 


IQ-22 


x 


Jh_ 

£  +  £i 


3 


0  <  £  <  0.59 Ih  —  £i 
£  >  0.59 Ih  —  £i 


(26) 


2.6  Bound-bound  transitions 


Bound-bound  transitions  are  another  significant  source  of  radiative  cooling  of  the  plasma. 
The  photo-absorption  cross  section  for  a  given  bound-bound  transition  may  be  estimated 
by  (Zel’dovich  &  Raizer), 


&ij  — 


au>ij  dv  =  2.65  x  10 


f? 

-6  J  l 


abs 

ij 


' iJ 


' iJ 


I  in 


(27) 


where  the  absorption  oscillator  strength  is  given  by 


[’abs 

ij 


ft  37  ' 


(28) 


Contributions  to  the  parameter  7  have  been  assumed  to  result  from  a  combination  of 
natural  and  pressure  line  broadening, 


T  Tnat  Tcol 


(29) 


given  by 


7nat 


8ir2e2v2 
3  me3 


=  2.47  x  10-2V  [s"1! 


(30) 


and 


7coi  =  —  =  2cr  vn 
Tco\ 


(31) 


respectively.  Based  on  the  resulting  mean-free  paths,  it  has  been  assumed  that  all  bound- 
bound  radiation  to  the  ground  level  from  the  excited  states  is  absorbed.  For  all  other 
transitions,  the  mean-free  path  is  several  orders  of  magnitude  greater  than  the  dimensions 
of  the  shock  tube  and  the  associated  radiation  is  assumed  to  escape. 


2.7  Bremsstrahlung  emission 


Free-free  transitions  have  been  incorporated  in  the  model  via  Kramer’s  formula  [37]  for 
Bremsstrahlung  emission.  Thus  the  corresponding  term  in  the  vector  of  source  terms  is: 


dEe  1  167T2  veZ2se6g 

dt  J  CR  3\/3  meh(47T£0c)3 

=  -1.42  x  10-40Ze2ffTe1/2n+ne  [J/m3  •  s] 


(32) 
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where  g  is  the  gaunt  factor  (here  taken  to  be  unity)  and  the  effective  charge  z*  * 
taken  to  be  1.67  to  better  match  the  experimental  conditions  given  by  Glass  &  Liu  [12]. 
Bremsstrahlung  emission  resulting  from  neutral  atoms  is  1-2  orders  of  magnitude  less  than 
for  ions  and  has  therefore  been  neglected  [37] .  The  plasma  is  assumed  to  be  optically  thin 
to  all  Bremsstrahlung  emission. 


2.8  Model  reduction 

Since  computational  costs  increase  exponentially  with  the  addition  of  more  states,  solving 
fully  detailed  CR  models  quickly  becomes  intractable,  especially  when  coupled  with  multi¬ 
dimensional  CFD  codes.  For  this  reason,  it  is  necessary  to  develop  reduced  CR  models 
from  which  accurate  kinetics  can  still  be  recovered  for  the  conditions  of  interest.  The 
reduced  models  must  be  benchmarked  against  the  most  detailed  models  available  to  ensure 
their  validity. 


3  Argon 

Due  to  its  relative  abundance,  argon  has  been  extensively  used  in  experimental  investi¬ 
gations  of  plasma  gas  dynamics.  Shocktube  experiments  in  particular  have  been  used  to 
validate  CR  models  (c.f.  [24,  23,  20,  27,  14]).  But  perhaps  the  most  significant  work  was 
carried  out  at  the  Institute  of  Aerospace  Sciences  at  the  University  of  Toronto  (UTIAS) 
in  which  unsteady  effects  were  discovered  to  be  a  further  characteristic  of  hypervelocity 
shocks  in  argon  (c.f.  [2,  13,  12]).  This  was  taken  advantage  of  by  Cambier  [5]  who  showed 
that  galloping  instabilities  of  the  shock  front  could  in  fact  be  explained  through  the  influ¬ 
ence  of  the  kinetics  of  excitation  and  ionization  on  the  transport  of  the  gas.  Subsequent 
work  by  Kapper  &  Cambier  [17,  18]  extended  the  sophistication  of  the  CR  model  and  were 
able  to  reproduce  the  unsteady  shock  structure  in  two-dimensional  space.  That  model  is 
reiterated  here  and  serves  as  the  foundation  for  extension  to  krypton  and  xenon. 

3.0.1  Atom  impact  processes 

For  atom-impact  excitation  of  argon,  we  have  used  Drawin’s  formula  (Eq.  (8))  for  the 
allowed  transitions.  The  inner  ‘ip' As  manifold  transitions,  however,  take  the  form  [4] 

<^)=Ph£-^L  (33) 

bij[eV] 

where  £ij[ev]  is  the  transition  energy  in  eV.  The  linear  /3-parameters  are  provided  in 
Tables  3.0.1  and  3.0.1.  Note  the  relative  strength  of  A jc  =  0  transitions  compared  to  the 
4s-4s'  (Ajcne0)  transition. 

For  atom-impact  ionization,  the  cross  section  for  ionization  from  ground  state  argon  is 
taken  from  Haugsjaa  and  Amme  [15], 

<r?(e)  =  1.8  x  10~25(£[eV]  -  15.760)1'3  [m2],  (34) 

while  the  formula  of  Drawin  (Eq.  10)  has  been  used  for  the  excited  states. 
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i 

e(i)  [eV] 

9i 

jc 

ni\K\j 

i 

e(i)  [eV] 

9i 

jc 

ni\K\j 

i 

0 

1 

1.5 

[Mg]3p6 

42 

14.711 

3 

1.5 

4d[l/2]i 

2 

11.548 

5 

1.5 

4s[3/2]2 

43 

14.738 

1 

0.5 

5p/[l/2]o 

3 

11.624 

3 

1.5 

4s[3/2]i 

44 

14.743 

5 

1.5 

4d[3/2]2 

4 

11.723 

1 

0.5 

4s'[l/2]0 

45 

14.757 

9 

1.5 

4d[7/2]4 

5 

11.828 

3 

0.5 

4s'[l/2]! 

46 

14.781 

7 

1.5 

4d[7/2]3 

6 

12.907 

3 

1.5 

4p[l/2h 

47 

14.809 

5 

1.5 

4d[5/2]2 

7 

13.076 

7 

1.5 

4p[5/2]3 

48 

14.824 

7 

1.5 

4d[5/2]3 

8 

13.095 

5 

1.5 

4p[5/2]2 

49 

14.839 

5 

1.5 

6s[3/2]2 

9 

13.153 

3 

1.5 

4p[3/2]i 

50 

14.848 

3 

1.5 

6s  [3/2] ! 

10 

13.172 

5 

1.5 

4p[3/2]2 

51 

14.859 

3 

1.5 

4d[3/2]4 

11 

13.273 

1 

1.5 

4p[l/2]o 

52 

14.901 

3 

1.5 

4/[3/2]4 

12 

13.283 

3 

0.5 

4p'[3/2]1 

53 

14.901 

5 

1.5 

4/ [3/2]  2 

13 

13.302 

5 

0.5 

4^  [3/2]  2 

54 

14.904 

11 

1.5 

4/ [9/2]  5 

14 

13.328 

3 

0.5 

4p'[l/2]1 

55 

14.904 

9 

1.5 

4/[9/2]4 

15 

13.480 

1 

0.5 

4p'[l/2]0 

56 

14.907 

7 

1.5 

4/ [5/2]  3 

16 

13.845 

1 

1.5 

3d[l/2]0 

57 

14.907 

5 

1.5 

4/ [5/2]  2 

17 

13.864 

3 

1.5 

3d[l/2]x 

58 

14.909 

7 

1.5 

4/[7/2]3 

18 

13.903 

5 

1.5 

3c7[3/2]2 

59 

14.909 

9 

1.5 

4/[7/2]4 

19 

13.979 

9 

1.5 

3d[7/2]4 

60 

14.953 

5 

0.5 

4d'  [3/2]  2 

20 

14.013 

7 

1.5 

3rf[7/2]s 

61 

14.955 

5 

0.5 

4d'[5/2]2 

21 

14.063 

5 

1.5 

3d[5/2]2 

62 

14.972 

7 

0.5 

4d'[5/2]3 

22 

14.068 

5 

1.5 

5s[3/2]2 

63 

15.004 

3 

0.5 

4d'[3/2]x 

23 

14.090 

3 

1.5 

5s[3/2]4 

64 

15.011 

3 

1.5 

6p[l/2]i 

24 

14.099 

7 

1.5 

3d  [5/2]  3 

65 

15.014 

1 

0.5 

6s'[l/2]0 

25 

14.153 

3 

1.5 

3d[3/2]i 

66 

15.022 

3 

0.5 

6s'  [1/2]  4 

26 

14.214 

5 

0.5 

3d' [5/2]  2 

67 

15.023 

7 

1.5 

6p[5/2]3 

27 

14.234 

5 

0.5 

3d' [3/2]  2 

68 

15.026 

5 

1.5 

6p[5/2]2 

28 

14.236 

7 

0.5 

3d' [5/2]  3 

69 

15.034 

3 

1.5 

6p[3/2]x 

29 

14.241 

1 

0.5 

5s'[l/2]0 

70 

15.036 

5 

1.5 

6p[3/2]2 

30 

14.255 

3 

0.5 

5s'[l/2]! 

71 

15.060 

1 

1.5 

6p[l/2]o 

31 

14.304 

3 

0.5 

3d' [3/2]  i 

72 

15.083 

7 

0.5 

4/' [7/2]  3 

32 

14.464 

3 

1.5 

5p[l/2]i 

73 

15.083 

9 

0.5 

4/'[7/2]4 

33 

14.499 

7 

1.5 

5p[5/2]3 

74 

15.083 

7 

0.5 

4/'[5/2]3 

34 

14.506 

5 

1.5 

5p[5/2]2 

75 

15.083 

5 

0.5 

4/' [5/2]  2 

35 

14.525 

3 

1.5 

5p[3/2]i 

76 

15.101 

1 

1.5 

5d[l/2]0 

36 

14.529 

5 

1.5 

5p[3/2]2 

77 

15.118 

3 

1.5 

5d[l/2]4 

37 

14.576 

1 

1.5 

5p[l/2]0 

78 

15.131 

9 

1.5 

5d[7/2]4 

38 

14.681 

3 

0.5 

5p'[3/2]! 

79 

15.137 

5 

1.5 

5d[3/2]2 

39 

14.687 

3 

0.5 

5p'[l/2]! 

80 

15.146 

7 

1.5 

5d[7/2]3 

40 

14.688 

5 

0.5 

5p'[3/2]2 

oo 

15.760 

4 

1.5 

[Mg]3p5 

41 

14.694 

1 

1.5 

4d[l/2]0 

oo' 

15.937 

2 

0.5 

[Mg]3p5 

Table  3:  Lowest  80  levels  of  Ar  I  by  energy. 
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i 

3 

3*. 

iij 

m 

2/eV] 

fij 

3 p6 

3d' [3/2]  i 

9.24 

X 

10-24 

1.06 

x  lO”01 

3  p6 

3d[3/2]i 

8.38 

X 

10-24 

9.32 

x  10_ 02 

3  p6 

4s'[l/2]1 

3.85 

X 

IQ-23 

2.50 

x  lO”01 

3  p6 

4s  [3/2]  x 

9.89 

X 

10-24 

6.09 

x  lO”02 

3 p6 

^[1/21 

1.05 

X 

10-24 

1.19 

x  10_ 02 

3 pe 

5s  [3/2] ! 

2.46 

X 

10-24 

2.70 

x  10_ 02 

4s'[1/2]0 

4p'[l/2]i 

5.38 

X 

IQ-22 

3.14 

x  lO”01 

4s'[l/2]0 

4p'[3/2]! 

9.88 

X 

10-22 

5.29 

x  lO”01 

4s'[l/2]0 

4p[l/2]i 

2.05 

X 

IQ-22 

4.80 

x  lO”02 

4s/[l/2]0 

4p[3/2]i 

1.99 

X 

IQ-22 

8.22 

x  lO"02 

4s/[l/2]1 

4p'[l/2]0 

1.96 

X 

10-22 

1.25 

x  lO”01 

4s'[l/2]1 

4p'[l/2]1 

3.30 

X 

10-22 

1.57 

x  lO”01 

4s'[l/2]1 

4p'[3/2]! 

3.48 

X 

IQ-22 

1.51 

x  lO”01 

4s'[l/2]1 

4^  [3/2]  2 

8.71 

X 

IQ-22 

3.94 

x  lO”01 

4s'[l/2]1 

4p[l/2]i 

2.14 

X 

IQ-23 

3.80 

x  lO"03 

4s/[l/2]1 

4p[3/2]i 

4.23 

X 

IQ-23 

1.39 

x  10_ 02 

4s/[l/2]1 

4p[3/2]2 

3.12 

X 

IQ-22 

1.07 

x  lO”01 

4s'[l/2]1 

4p[5/2]2 

1.23 

X 

IQ-22 

3.52 

x  lO”02 

4s  [3/2]  x 

4p'[l/2]0 

5.82 

X 

IQ-25 

5.26 

x  lO"04 

4s  [3/2]  x 

4p'[l/2]! 

2.07 

X 

IQ-23 

1.45 

x  10_ 02 

4s  [3/2]  x 

4p/[3/2]1 

2.79 

X 

IQ-25 

1.80 

X  10“°4 

4s  [3/2] ! 

4p'[3/2]2 

1.72 

X 

IQ-22 

1.15 

x  lO”01 

4s  [3/2]  x 

4p[l/2]0 

1.80 

X 

IQ-22 

1.14 

x  lO”01 

4s  [3/2]  x 

4p[l/2]i 

2.55 

X 

IQ-22 

7.60 

X  10-°2 

4s  [3/2]  x 

4p[3/2]i 

4.95 

X 

IQ-22 

2.50 

x  lO”01 

4s  [3/2] ! 

4p[3/2]2 

1.50 

X 

IQ-22 

7.85 

x  10_ 02 

4s  [3/2] ! 

4p[5/2]2 

8.47 

X 

IQ-22 

3.81 

x  lO”01 

4s  [3/2]  2 

4p'[l/2]i 

3.51 

X 

IQ-23 

2.79 

x  10_ 02 

4s  [3/2]  2 

4p/  [3/2]: 

3.90 

X 

lO-24 

2.87 

x  lO”03 

4s  [3/2]  2 

4^  [3/2]  2 

3.74 

X 

IQ-23 

2.85 

x  lO”02 

4s  [3/2]  2 

4p[l/2]i 

4.01 

X 

IQ-22 

1.42 

x  lO”01 

4s  [3/2]  2 

4p[3/2]i 

4.76 

X 

IQ-23 

2.78 

x  10_ 02 

4s  [3/2]  2 

4p[3/2]2 

3.54 

X 

IQ-22 

2.14 

x  lO”01 

4s  [3/2]  2 

4p[5/2]2 

1.71 

X 

IQ-22 

8.94 

x  lO”02 

4s  [3/2]  2 

4p[5/2]3 

9.10 

X 

IQ-22 

4.58 

x  lO”01 

Table  4:  Atom  impact  excitation  parameters  for  allowed  transitions  from  ground  state 
Ar. 
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Figure  1:  Partial  Grotrian  diagram  for  neutral  argon  taken  from  [30]. 


i 

j 

ftj  [m2/eV] 

4s  [3/2]  2 

4s  [3/2] ! 

1.79  x  10"24 

4s  [3/2]  2 

4s,[l/2]0 

4.80  x  10-26 

4s  [3/2]  2 

4s/[l/2]1 

4.80  x  10-26 

4s  [3/2] ! 

4s'[l/2]0 

4.80  x  10"26 

4s  [3/2] ! 

4s'[l/2]1 

4.80  x  10"26 

4s'[l/2]0 

4s'[l/2]1 

1.79  x  10"24 

Table  5:  Atom  impact  excitation  parameters  for  intra-4s  transitions. 
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3.1  Electron  impact  processes 

Zatsarinny  &  Bartschat  [35]  have  performed  A-matrix  calculations  to  obtain  highly  de¬ 
tailed  cross  sections  for  electron  impact  excitation  of  argon  near  threshold  energies.  These 
have  been  provided  for  all  inter-level  transitions  up  to  the  3d  and  5s  manifolds.  Beyond 
this,  Drawin’s  formula  has  been  used  systematically. 

For  ionization  from  the  excited  states,  we  have  utilized  the  Deutsch-Mark  formalism  as 
previously  mentioned.  The  necessary  parameters  for  Eqs.  (12)  and  (13)  are  given  in 
Table  3.1,  while  the  resulting  cross  sections  have  been  plotted  in  Figure  2.  For  ionization 
from  the  excited  states  we  have  considered  the  electron  occupying  the  valence  shell  only. 
For  ionization  from  the  core,  we  have  considered  contributions  form  all  sub-shells.  For 
ionization  from  the  core,  we  have  have  compared  the  cross  sections  from  the  DM  formalism 
with  the  experimental  fit  of  Straub  et  al  [32] . 


3.2  Elastic  collisions 

The  theoretical  cross  section  of  McEachran  &  Stauffer  [21]  have  been  utilized  for  the 
elastic  momentum  transfer  cross  sections  between  neutral  argon  and  electrons,  which  are 
reproduced  in  Table  3.2  and  plotted  in  Figure  3. 


3.3  Photorecombination/Photoionization 

For  argon,  the  parameter  7 %  in  Eq.  (26)  takes  the  values  0.0763,  0.0458,  0.0305,  and  0.0915 
for  i  =  2,  3,  4,  and  5,  respectively  [33].  Based  on  nominal  shock  tube  plasma  conditions, 
an  estimation  of  the  mean-free  path  for  photoionization  from  the  ground  state  is  on  the 
order  of  100  pm.  In  the  case  of  the  4s  levels,  their  lower  populations  combined  with 
reduced  cross  sections  result  in  mean-free  paths  on  the  order  of  1  m.  Consequently,  free- 
bound  radiation  into  the  ground  state  is  assumed  to  be  locally  absorbed,  while  the  plasma 
is  assumed  to  be  optically  thin  to  transitions  from  the  excited  levels. 


Valence  electron 
shell,  n£ 

T'ni  [A] 

9ni  X  In£  [gV] 

4  s 

2.49 

7.40 

5  s 

6.35 

6.35 

4 p 

3.40 

31.00 

3d 

4.36 

13.60 

Table  6:  Radii  of  Ar  valence  electron  and  reduced  weighting  factors  for  £  =  1  as  taken 
from  [8]  and  [6]. 
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Figure  2:  Ionization  cross  sections  due  to  electron  impact  for  Ar  I  as  computed  by  the 
Deutsch-Mark  formalism. 


e  [eV] 

^  [1020  m2] 

e  [eV] 

den  [1020  m2] 

0.01 

4.0231 

0.41 

0.3396 

0.03 

2.3425 

0.51 

0.5349 

0.05 

1.5247 

0.61 

0.7381 

0.07 

1.0261 

0.71 

0.9318 

0.09 

0.6998 

0.81 

1.1139 

0.13 

0.3291 

0.91 

1.2835 

0.17 

0.1599 

1.00 

1.4280 

0.19 

0.1198 

1.50 

2.1438 

0.21 

0.0983 

2.00 

2.8318 

0.23 

0.0917 

3.00 

4.4046 

0.25 

0.0969 

4.00 

6.3502 

0.29 

0.1316 

5.00 

8.6317 

0.31 

0.1573 

7.50 

14.6559 

0.32 

0.1724 

10.00 

17.8325 

Table  7:  Momentun  transfer  cross  sections  for  electron-neutral  collisions  [21]. 


Figure  3:  Momentum  transfer  cross  sections  for  Ar  I. 
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lower 

level 

upper 

level 

mean-free 
path  [m] 

3 p6 

4s  [3/2] ! 

3.59  x  10-2 

3 p6 

Ml/2]! 

9.47  x  10"3 

3 p6 

3d[l/2]x 

2.41 

3 p6 

5s  [3/2] ! 

7.11  x  10-2 

3 p6 

3d  [3/2]  i 

2.54  xl0“2 

3 p6 

5VI1/2]! 

1.34  xKT1 

3 p6 

3d' [3/2]! 

3.05  x  10"2 

Table  8:  Estimated  mean- free  paths  for  bound-bound  transitions  to  ground  state  for 
T  =  1  eV  and  n3pe  =  1024  m-3. 

3.4  Bound-bound  transitions 


Bound-bound  transitions  are  another  significant  source  of  radiative  cooling  of  the  plasma. 
The  photo-absorption  cross  section  for  a  given  bound-bound  transition  may  be  estimated 
by  [37], 


1  7  f- 

<jij  =  —  /  <jtj_a  du  =  2.65  x  10-6— 


abs 


V, 


m 


(35) 


where  the  absorption  oscillator  strength  is  given  by 


Tabs 

J  ij 


9j  Ajj 

9i  37  ’ 


(36) 


For  argon,  we  have  estimated  the  contributions  to  the  parameter  7  in  Eq.  28  by  assuming 
a  combination  of  natural  and  pressure  line  broadening, 


T  Tnat  "F  Tcol 


(37) 


given  by 


and 


7nat 


8n2e2v2 

8m<F 


=  2.47  x  10-2V  is"1] 


7coi  =  —  =  2  avn 
Tcol 


(38) 

(39) 


respectively.  The  collision  time  rcoi  and  the  collision  cross  section  a  have  been  tallied  for 
all  argon-argon  collisions.  The  radiation  cross  sections  Eq.  (35)  have  been  estimated  using 
the  oscillator  strengths  and  transition  probabilities  from  Ralchenko  et  al  [29]  along  with 
representative  plasma  conditions  of  the  shock  tube  experiments.  The  resulting  mean- free 
paths  for  bound-bound  transitions  to  ground  state  given  in  Table  VII  are  below  the  length 
scales  of  the  shock  tube  dimensions  and  induction  length  for  most  of  the  levels. 
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3.5  Validation 

This  current  model  has  been  extensively  benchmarked  in  the  case  of  shock-heated  argon 
plasma.  Results  obtained  by  Kapper  &  Cambier  are  shown  in  Figure  4.  A  detailed 
analysis  of  the  results  are  provided  within. 


4  Krypton 

Beyond  argon,  the  heavier  noble  gases  including  krypton  and  xenon  present  a  considerable 
challenge  to  theoretical  computation  of  oscillator  strengths  and  subsequent  particle  impact 
cross  sections.  The  significant  influence  of  relativistic  effects  means  that  models  developed 
for  lighter  species  have  been  unsuccessful  in  obtaining  quantitative  agreement.  Despite 
this,  active  progress  is  being  made  in  extending  the  models  that  will  greatly  benefit  our 
CR  models. 

The  lowest  80  levels  of  Kr  I  as  given  by  NIST  are  listed  in  Table  9.  Note  that  the  energy 
needed  to  excite  an  electron  to  the  lowest  manifold  is  approximately  1.5  eV  lower  than 
that  or  Ar  I,  with  subsequent  ionization  requiring  an  additional  4  eV,  similar  to  argon. 
This  translates  to  higher  ionization  fractions  behind  the  shock. 

4.0.1  Atom  impact  processes 

For  excitation  due  to  heavy  particle  impact,  the  formula  of  Drawin  (Eq.  8)  has  been  used 
to  obtain  cross  sections  for  the  allowed  transitions.  The  parameters  B*-  are  provided  in 
Table  4.0.1. 


4.1  Electron  impact  processes 

Building  upon  their  success  for  electron-impact  excitation  cross  sections  for  argon,  [36] 
have  extended  their  calculations  to  krypton  and  xenon.  The  computed  cross  sections 
reproduced  in  Figures  10  and  11  show  good  agreement  with  high- resolution  experiments 
near  threshold  for  transitions  to  the  5s  manifold  from  the  ground  state. 

For  ionization  from  the  excited  states,  we  have  used  the  Deutsch-Mark  formalism  and 
have  considered  the  electron  in  the  valence  shell  only  since  this  leads  to  better  agreement 
with  experimental  data  [7].  For  ionization  from  the  core,  we  have  considered  only  the  5s 
sub-shell.  The  resulting  cross  sections  have  been  plotted  in  Figure  6. 


4.2  Elastic  collisions 

The  momentum  transfer  cross  sections  for  elastic  impact  of  electrons  with  krypton  have 
been  given  by  Mitroy  [25]  for  energies  up  to  1  eV.  Beyond  this,  we  have  used  the  cross 
sections  of  McEachran  &  Stauffer  [22],  The  resulting  cross  sections  are  plotted  in  Figure 
7. 
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Figure  4:  Heavy  particle  mass  density  and  electron  number  density  profiles  (top)  and 
ionization  fractions  ( bottom )  behind  Mach  15.9  (left)  and  Mach  16.1  (right)  shocks  in 
pure  argon.  Results  shown  have  been  obtained  considering  excited  levels  from  the  4s,  4 p, 
5s,  and  3d  manifolds. 


Figure  5:  Electron- impact  excitation  cross  sections  for  krypton  taken  from  [36]. 
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Table  9:  Lowest  80  levels  of  Kr  I  by  energy. 
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Table  10:  Atom  impact  excitation  parameters  for  allowed  transitions  from  ground  state 

Kr  I. 
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Figure  6:  Ionization  cross  sections  due  to  electron  impact  for  Kr  I  as  computed  by  the 
Deutsch-Mark  formalism. 


Figure  7:  Momentum  elastic  cross  sections  for  Kr  I. 

4.3  Line  transitions 

The  Einstein  coefficients  and  oscillator  strengths  for  allowed  transitions  available  through 
the  NIST  database  are  provided  in  Table  4.3.  The  two  transitions  from  ground  state 
krypton  have  been  computed  by  Wilkinson  [34] . 

4.4  Validation 

Initial  benchmarks  of  the  krypton  CR  model  have  been  carried  out  for  the  experimental 
conditions  provided  by  Glass  et  al  [13].  Preliminary  results  are  provided  in  Figure  8  for 
which  only  excited  levels  of  the  5s  manifold  have  been  considered.  Furthermore,  radiation 
effects  have  been  neglected. 


5  Xenon 

Xenon  is  the  preferred  propellant  in  ion  thrusters  including  Hall  thrusters  due  to  its  high 
Z- value  and  resulting  /sp.  Due  to  the  highly  nonequilibrium  conditions  under  which  such 
devices  operate,  a  detailed  CR  model  is  necessary  to  model  the  flow  dynamics.  Computing 
the  required  cross  sections,  however,  is  a  challenge  due  to  relativistic  effects  as  mentioned 
for  krypton.  As  such,  the  models  that  have  been  successfully  applied  to  argon  may 
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Table  11:  Allowed  transitions  for  Kr  I  available  from  NIST  [29].  The  oscillator  strength 
is  denoted  by  /. 
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Figure  8:  Heavy  particle  mass  density  and  electron  number  density  profiles  (top)  and 
ionization  fractions  ( bottom )  behind  Mach  15.17  (left)  and  Mach  15.05  (right)  shocks  in 
pure  krypton.  Results  shown  have  been  obtained  considering  only  excited  levels  from  the 
5s  manifold  while  neglecting  radiative  processes. 


breakdown  and  therefore  caution  must  be  exercised  in  their  application.  These  models 
have  been  used  as  an  initial  base  to  construct  out  CR  model  and  new  quantum  calculations 
which  do  take  into  account  the  necessary  relativistic  effects  are  promising.  Our  preliminary 
model  for  Xenon  follows. 

The  effect  of  spin-orbit  splitting  of  the  ion  due  to  the  core  angular  momentum  jc  (3/2  or 
1/2)  is  more  evident  in  Xe  than  Ar  which  leads  to  larger  energy  separation  between  the 
two  groups.  This  effect  can  be  seen  upon  comparison  of  the  Grotrian  diagrams  for  Xe 
and  Ar  in  Figure  9.  The  first  80  levels  are  provided  in  Table  12. 


5.1  Atom  impact  processes 

As  for  argon  and  krypton,  the  formula  of  Drawin  (Eq.  8)  has  been  used  to  obtain  cross 
sections  for  the  allowed  transitions.  The  parameters  /3f  are  provided  in  Table  5.1. 


5.2  Electron  impact  processes 

Zatsarinny  &  Bartschat  [36]  have  computed  preliminary  electron  impact  excitation  cross 
sections  for  xenon.  Cross  sections  taken  from  their  publications  reproduced  in  Figures  10 
and  11  show  good  agreement  with  high-resolution  experiments  near  threshold. 
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Figure  9:  Comparison  of  Xe  and  Ar  energy  lines  from  [16]. 


Figure  10:  e- impact  excitation  cross  sections  to  5s  manifold  from  ground  state  Xe. 
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1.5 

8s  [3/2]! 

7 

9.686 

5 

1.5 

H  5/2]2 

48 

11.276 

7 

1.5 

4/[7/2]3 

8 

9.721 

7 

1.5 

H  5/2]3 

49 

11.276 

9 

1.5 

4/[7/2]4 

9 

9.789 

3 

1.5 

6p[3/2]i 

50 

11.301 

5 

0.5 

5(f  [5/2]2 

10 

9.821 

5 

1.5 

6p[3/2]2 

51 

11.338 

5 

0.5 

5(f  [3/2]2 

11 

9.890 

1 

1.5 

5d[l/2]0 

52 

11.375 

7 

0.5 

5(f  [5/2]3 

12 

9.917 

3 

1.5 

5rf[l/2]x 

53 

11.422 

3 

1.5 

7d[l/2]i 

13 

9.933 

1 

1.5 

6p[l/2]0 

54 

11.426 

3 

1.5 

8p[l/2]i 

14 

9.943 

9 

1.5 

5d[7/2]4 

55 

11.434 

5 

1.5 

8p[5/2]2 

15 

9.959 

5 

1.5 

5d[3/2]2 

56 

11.439 

1 

1.5 

7d[l/2]0 

16 

10.039 

7 

1.5 

5d[7/2]3 

57 

11.439 

7 

1.5 

8p[5/2]3 

17 

10.157 

5 

1.5 

5d[5/2]2 

58 

11.448 

3 

1.5 

8p[3/2]i 

18 

10.220 

7 

1.5 

5d[5/2]3 

59 

11.453 

5 

1.5 

8p[3/2}2 

19 

10.401 

3 

1.5 

5rf[3/2]i 

60 

11.462 

9 

1.5 

7d[7/2}4 

20 

10.562 

5 

1.5 

7s[3/2]2 

61 

11.475 

1 

1.5 

8p[l/2]0 

21 

10.593 

3 

1.5 

7s  [3/2]! 

62 

11.487 

7 

1.5 

7d[7/2]3 

22 

10.902 

3 

1.5 

7p[l/2]i 

63 

11.491 

5 

1.5 

7c?[5/2]2 

23 

10.954 

5 

1.5 

7p[5/2  ]2 

64 

11.495 

3 

1.5 

7d[3/2]i 

24 

10.958 

3 

0.5 

6p'[3/2]! 

65 

11.496 

5 

1.5 

7d[3/2]2 

25 

10.969 

7 

1.5 

7p[5/2]3 

66 

11.498 

7 

1.5 

7d[5/2]3 

26 

10.971 

1 

1.5 

6d[l/2]0 

67 

11.575 

3 

1.5 

5/[3/2]i 

27 

10.979 

3 

1.5 

6d[l/2]x 

68 

11.576 

5 

1.5 

5/ [3/2]  2 

28 

10.996 

5 

1.5 

7p[3/2]2 

69 

11.577 

11 

1.5 

5/ [9/2]  5 

29 

10.998 

5 

1.5 

6d[3/2]2 

70 

11.577 

9 

1.5 

5/[9/2]4 

30 

11.003 

3 

1.5 

7p[3/2]i 

71 

11.580 

5 

1.5 

9s[3/2]2 

31 

11.015 

1 

1.5 

7p[l/2]0 

72 

11.580 

7 

1.5 

6/ [5/2]  3 

32 

11.024 

9 

1.5 

6d[7/2]4 

73 

11.581 

5 

1.5 

5/ [5/2]  2 

33 

11.038 

7 

1.5 

6d[7/2]3 

74 

11.583 

7 

1.5 

5/[7/2]3 

34 

11.055 

5 

0.5 

6p'[3/2]2 

75 

11.583 

9 

1.5 

5/[7/2]4 

35 

11.065 

5 

1.5 

6d[5/2]2 

76 

11.583 

5 

1.5 

5^  [5/2]  2 

36 

11.069 

3 

0.5 

6p,[l/2]1 

77 

11.583 

7 

1.5 

5s1  [5/2]  3 

37 

11.101 

7 

1.5 

6d[5/2]3 

78 

11.583 

3 

1.5 

9s  [3/2]  i 

38 

11.141 

1 

0.5 

6p'[l/2]0 

79 

11.584 

11 

1.5 

5^[ll/2]5 

39 

11.163 

3 

1.5 

6d[3/2]1 

80 

11.584 

13 

1.5 

5fl,[H/2]6 

40 

11.258 

5 

1.5 

8s[3/2]2 

oo 

12.130 

4 

1.5 

[Cd]5p5 

41 

11.263 

3 

1.5 

4/[3/2]x 

oo' 

13.436 

2 

0.5 

[Cd]5 p5 

Table  12:  Lowest  80  levels  of  Xe  I  by  energy. 
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Figure  11:  e-impact  excitation  cross  sections  to  5 p  manifold  from  ground  state  Xe  ([1]). 
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i 

j 

d* 

'  *7 

m2/eV] 

fij 

5  p6 

5d'[3/2]i 

3.34 

x  10-23 

2.05 

x  10“01 

5 p6 

5d[l/2]x 

2.56 

x  10-24 

9.80 

x  10“03 

5 p6 

5d[3/2]i 

9.06 

x  10“23 

4.00 

x  10“01 

5  p6 

6d[l/2]i 

4.82 

x  10-25 

2.50 

x  10“03 

5  p6 

6d[3/2]! 

1.52 

x  10-23 

8.30 

x  10-°2 

5  p6 

Ss'll/2]! 

5.41 

x  10-23 

1.86 

x  10“01 

5 p6 

6s  [3/2]  i 

1.16 

x  10-22 

2.73 

x  10“01 

5 p6 

7d[l/2]x 

3.71 

x  10-24 

2.17 

x  10-°2 

5 p6 

7s  [3/2]  x 

2.12 

x  10“23 

9.90 

x  10-°2 

5  p6 

8s[3/2]x 

4.20 

x  10-24 

2.36 

x  10-°2 

6p[  5/2]3 

7d[7/2]4 

8.73 

x  10-23 

6.50 

x  10-°2 

6s'[l/2]0 

7p[3/2]x 

1.62 

x  10-23 

8.60 

x  10-°3 

6s[3/2]x 

7p[l/2]0 

1.16 

x  10-24 

2.80 

x  10“03 

6s  [3/2]  2 

6p'[l/2]x 

3.73 

x  10“25 

1.10 

x  10“03 

6s  [3/2]  2 

6^  [3/2]  2 

2.21 

x  10-25 

6.40 

x  10“04 

6s  [3/2]  2 

6p[  3/2] ! 

1.42 

x  10-23 

6.40 

x  10“03 

6s  [3/2]  2 

7p[3/2]2 

8.46 

x  10-25 

2.30 

x  10“03 

6s  [3/2]  2 

7p[5/2]3 

1.74 

x  10-24 

4.60 

x  10“03 

Table  13:  Atom  impact  excitation  parameters  for  allowed  transitions  of  Xe  I. 


5.3  Elastic  collisions 


The  momentum  transfer  cross  sections  for  Xe  have  been  plotted  in  Figure  12  [26].  These 
data  include  cross  sections  from  [31]  for  energies  up  to  2  eV  and  from  Gibson  et  al  [11] 
for  energies  up  to  50  eV.  The  values  of  Gibson  were  extrapolated  to  obtain  cross  sections 
to  100  eV. 


Figure  12:  Momentum  elastic  cross  sections  for  Xe  I  compiled  by  Mitroy  [26]. 
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5.4  Line  transitions 

Like  Krypton,  only  a  handful  of  allowed  transitions  for  Xe  I  are  provided  in  the  NIST 
database  which  are  given  in  Table  5.4.  As  previously  mentioned,  relativistic  effects  add 
complexity  to  ab  initio  calculations,  resulting  in  significant  uncertainty. 

5.5  Validation 

While  simulations  in  Xenon  have  yet  to  be  performed,  we  present  here  experimental  data 
for  future  benchmarks  of  the  model.  The  experimental  data  obtained  by  Ezumi  et  al  [9] 
for  a  Mach  13.1  shock  in  pure  xenon  is  plotted  in  Figure  13.  Upon  fitting  his  model  to 
the  data,  Ezumi  was  able  to  determine  a  value  of  1.8  x  10-13  m2/eV  for  the  xenon-xenon 
excitation  cross  section.  Data  has  also  been  provided  for  hydrogen-xenon  mixtures,  from 
which  the  corresponding  atom-atom  excitation  cross  sections  can  be 


6  Rate  formulation 

The  cross  sections  must  be  integrated  based  on  the  distribution  of  particle  energies  in  order 
to  be  applied  successfully  in  a  meaningful  computation.  The  forward  rate  coefficients  are 
computed  assuming  a  Maxwellian  electron  energy  distribution  function  (EEDF)  according 
to 


v  r°° 

k  =  (ksT) 2  J  a^£e~£/kBT  d£  (40) 

where  e0  is  the  threshold  energy  and  v  is  the  mean  thermal  velocity, 

8  fcBT\1/2 

7T^  J 

with  n  the  reduced  mass.  The  backward  rates  are  then  computed  from  the  principle  of 
detailed  balance.  For  excitation  processes  this  takes  the  form 


(41) 


h,dex  _  I  exc9i  ea/kBT 

3  13  9j 

where  £%j  =  £j  —  e-L  is  the  energy  difference  between  the  upper  and  lower  levels.  The 
temperature  in  Eq.  (42)  is  that  of  the  the  impact  species.  For  ionization  and  recombination 
processes,  the  reverse  rates  are  also  computed  via  detailed  balance,  leading  to: 


Zi 

urec  _  uton  1 

ki  ~ki  z+z„ 


where  Ze  is  the  electron  partition  function 


(43) 
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line  [nm 

A  [s  1 

/ 

lower  level 

upper  level 

104.38 

5.90xl0+°7 

2.90xl0-02 

5 p6 

7s/[l/2]1 

104.71 

1.30  x  10+08 

6.40xl0-02 

5 p6 

9d[3/2]i 

105.01 

8.50  xlO+06 

4.20xl0-03 

5 p6 

9d[l/2]i 

105.61 

2.45  xlO+08 

1.23xlO-01 

5 p6 

8d[3/2]1 

106.13 

1.90xl0+°7 

9.60  xlO"03 

5 p6 

8d[  1/2]! 

106.82 

3.99  xlO+08 

2.05  xlO“01 

5 p6 

5d,[3/2]1 

108.54 

4.10xl0+°7 

2.17xlO-02 

5 pe 

7d[l/2]i 

109.97 

4.34xlO+07 

2.36xlO-02 

5 p6 

8s  [3/2] ! 

111.07 

1.50  x  10+08 

8.30xl0-02 

5 p6 

6d[3/2]i 

112.93 

4.40  xlO+06 

2.50xl0-03 

5 p6 

6d[l/2]! 

117.04 

1.60xl0+°8 

9.90xl0-02 

5 p6 

7s[3/2]i 

119.20 

6.20  xlO+08 

4,OOxlO-01 

5 pe 

5d[3/2]i 

125.02 

1.40  xlO+07 

9.80  xlO"03 

5 p6 

5d[l/2]i 

129.56 

2.46  xlO+08 

1.86xlO-01 

5 p6 

6s,[l/2]1 

146.96 

2.81  xlO+08 

2.73xlO-01 

5 p6 

6s[3/2]i 

450.10 

6.20xl0+°5 

l.lOxlO-03 

6s  [3/2]  2 

6p'[l/2]! 

452.47 

2.10xl0+°5 

6.40  xlO“04 

6s  [3/2]  2 

6p'[3/2]2 

462.43 

7.20  xlO+05 

2.30xl0-03 

6s[3/2]2 

7p[3/2]2 

467.12 

1.00xl0+06 

4.60  xlO"03 

6s[3/2]2 

7p[5/2]3 

480.70 

2.40  xlO+06 

2.80xl0-03 

6s  [3/2] ! 

7p[l/2]0 

711.96 

6.60  xlO+06 

6.50xl0-02 

6p[5/2]3 

7d[7/2]4 

796.73 

3.00  xlO+05 

8.60  xlO“03 

6s'[l/2]0 

7p[3/2]i 

840.92 

1.00xl0+°6 

6.40  xlO“03 

6s[3/2]2 

6p[3/2]1 

Table  14:  Allowed  transitions  for  Xe  I  available  from  NIST  [29].  The  oscillator  strength 
is  denoted  by  /. 
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Figure  13:  Number  density  profiles  for  atoms  and  electrons  (top)  and  ionization  fractions 
(bottom)  behind  a  Mach  13.1  shock  (p  =  2.0  torr).  Left  profiles  are  for  pure  Xe,  while 
right  profiles  are  for  driven  gas  doped  with  0.17%  H2. 


(  27 imekBTe  \  3//2 
= 2  (-^-) 

Since  the  electronic  states  are  considered  separately,  the  ratio  of  partition  functions  for 
Ar  and  Ar+  is  simply 


~  9i  Ph/KT 

t;  e 

Z+  9+ 


(45) 


with  Ii  being  the  ionization  potential  of  the  ith  excited  state,  allowing  Eq.  (43)  to  be 
written  as 


urec  _  uion  9i  ^ 

'  9+  2 


K 2 


2'KmPknTP 


3/2 


Ji/KT 


(46) 


Note  that  the  exponential  temperature  dependence  in  Eq.  (46)  is  a  function  of  the  third 
body — T  is  taken  to  be  either  the  electron  or  heavy  particle  temperature  for  cases  in 
which  the  third  body  is  an  electron  or  atom,  respectively. 

In  addition  to  species  production  rates,  the  second  moment  with  respect  to  energy  is 
also  necessary  to  determine  the  energy  production  rate  for  photorecombination  (radiative 
capture)  processes, 
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k = v^v  r  a(-e'>e2e~‘,t‘T  *•  <47> 

In  practice,  all  rate  coefficients  and  their  derivatives  are  computed  a  priori  and  tabulated 
as  a  function  of  temperature. 

The  model  cross  sections  and  rates  used  in  the  CR  model  are  now  described.  For  further 
details,  the  reader  is  directed  to  the  work  of  Bultel  et  al  [4]  and  Vlcek  [33]  and  the 
references  therein. 


7  Database  implementation 


The  collisional-radiative  models  and  theory  of  the  preceding  sections  have  been  incorpo¬ 
rated  into  a  database  framework  that  can  manage  and  manipulate  cross  sections.  Dubbed 
“Colorado”  for  Collisional  Radiative,  the  database  serves  three  main  purposes.  The  first 
is  a  management  system  which  organizes  the  multitude  of  cross  sections  in  an  easily  ac¬ 
cessible  fashion.  The  second  is  to  provide  a  means  to  compute  cross  sections  based  on 
semi-empirical  formulas  described  in  the  preceding  sections  when  no  better  alternatives 
are  available.  The  third,  and  perhaps  most  critical  function,  is  to  compute  reaction  rates 
from  the  aggregate  of  cross  sections,  yielding  a  detailed  kinetics  database  that  is  directly 
accessible  to  CFD  solvers. 

Furthermore,  Colorado  is  an  object-orientated  (00)  framework  that  makes  use  of  00 
constructs  to  simplify  the  hierarchical  complexity  of  state-specific  nonequilibrium  within 
atoms  and  molecules.  Fine-grained  00  design  allows  encapsulation  of  necessary  data  for 
electronic  levels  and  transition  lines  for  a  more  physical  representation  of  the  data,  while 
coarse-graining  of  computational-intensive  tasks  preserves  efficiency. 


7.1  Structures 

From  a  user  standpoint,  the  code  is  fairly  straight-forward.  A  user  provides  an  input  file 
that  specifies  the  levels  for  which  the  kinetic  rates  should  be  computed.  Here,  however, 
we  go  a  step  beyond  what  is  necessary  for  execution  of  the  database  and  discuss  some  of 
the  inner  workings. 

When  the  database  is  executed,  information  about  the  electronic  levels  and  line  transitions 
are  read  in  from  existing  files  for  a  given  element.  These  data  form  a  collection  of  levels  and 
lines  from  which  the  user  can  extract  a  specific  subset,  based  on  a  simple  input  file.  The 
database  then  performs  tasks  including  cross  section  calculations  and  rate  computations 
for  the  selected  subset  of  levels  and  lines.  We  describe  the  treatment  for  the  levels  and 
lines  in  the  following  sections. 
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Configuration 

1 

1  . 

Term 

1 

-  1 

J  1 
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Level 

1 

1 

3s2 . 3p6 

1 

1 

1 

1 

IS 
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1 
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1 
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1 
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0.00 

1 

1 
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1 
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1 

1 

2  [3/2]* 

1 

1 

1 
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1 

1 

1 

1 

1 

1 
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1 
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1 

1 
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1 

1 
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1 

1 

1 

0  I 
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1 

1 

1 

1 

1 

1 
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i 
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1 

1 
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1 

1 

1 

2 [1/2] 

1 

1 

1 

1 
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1 
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1 

1 

1 
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1 

1 
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1 

1 

l 

3  I 
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1 

1 

1 

1 

1 

1 

2  I 

1 

13.09487173 

1 

1 

3s2 . 3p5 . (2P*<3/2>) .4p 

1 

1 
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1 

1 

1 

1  1 
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1 

1 

1 

1 

1 

1 

2  | 
l 

13.17177687 

1 

1 

3s2 . 3p5 . (2P*<3/2>) .4p 

1 

1 

1 

2 [1/2] 

1 

1 

1 

1 

0  I 

1 

13.27303726 

1 

1 

1 

3s2.3p5. (2P*<l/2>) .4p 

1 

1 

2 [3/2] 

1 

1 

1 

1  | 

13.28263817 

1 

1 

1 

1 

1 

1 

2  | 

1 

13.30222662 

1 

1 

3s2.3p5. (2P*<l/2>) .4p 

1 

1 

2 [1/2] 

1 

1 

1 

1  | 

13.32785620 

1 

1 

1 

1 

1 

1 

0  I 

i 

13.47988596 

1 

1 

3s2.3p5. (2P*<3/2>) .3d 

1 

1 
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1 

1 

1 

0  1 

13.8450376 

1 

1 

1 

1 

1 

1 

1  1 

1 

13.8636677 

1 

1 
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1 

1 

1 
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1 

1 

1 

1 

2  | 
l 
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1 

1 

1 
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1 

1 
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1 

1 

1 

4  | 

13.9792365 

1 

1 

1 

1 

1 

1 

3  | 

l 

14.0127372 

1 

1 

3s2 . 3p5 . (2P*<3/2>) . 3d 

1 

1 

1 

2  [5/2]* 

1 

1 

1 

1 

2  | 
l 

14.0630263 

1 

1 

1 

3s2 . 3p5 . (2P*<3/2>) . 5s 

1 

1 

2  [3/2]* 

1 

1 

1 

2  | 

14.0682968 

1 

1 

1 

1 

1 

1 

1  1 

1 

14.0899676 

1 

1 

3s2.3p5. (2P*<3/2>) .3d 

1 

1 

1 

2  [5/2]* 

1 

1 

1 

1 

3  I 

i 

14.0990550 

1 

1 

1 

3s2.3p5. (2P*<3/2>) .3d 

1 

1 

1 

2  [3/2]* 

1 

1 

1 

1 

1  | 
i 

14.1525142 

1 

1 

1 

3s2.3p5. (2P*<l/2>) .3d 

1 

1 

1 

2  [5/2]* 

1 

1 

1 

1 

2  | 

1 

14.2136706 

1 

1 

1 

3s2.3p5. (2P*<l/2>) .3d 

1 

1 

1 

2  [3/2]* 

1 

1 

1 

1 

2  | 

1 

14.2340217 

1 

1 

1 

3s2.3p5. (2P*<l/2>) .3d 

1 

1 

1 

2  [5/2]* 

1 

1 

1 

1 

3  I 

1 

14.2361052 

1 

1 

1 

3s2 . 3p5 . (2P*<l/2>) . 5s 

1 

1 

2  [1/2]* 

1 

1 

1 

0  I 

14.2410268 

1 

1 

1 

1 

1 

1 

1  1 

1 

14.2550847 

1 

1 

3s2.3p5. (2P*<l/2>) .3d 

1 

1 

1 

2  [3/2]* 

1 

1 

1 

1 

1  1 

1 

14.3036675 

1 

1 

1 

3s2 . 3p5 . (2P*<3/2>) . 5p 

1 

1 

1 

2 [1/2] 

1 

1 

1 

1 

1  1 
i 

14.46399485 

1 

1 

1 

3s2 . 3p5 . (2P*<3/2>) . 5p 

1 

1 

2  [5/2] 

1 

1 

1 

3  I 

14.49905272 

1 

1 

1 

1 

1 

1 

2  I 

i 

14.50606672 

1 

1 

3s2 . 3p5 . (2P*<3/2>) . 5p 

1 

1 

2 [3/2] 

1 

1 

1 

1  | 

14.52491238 

1 

1 

1 

1 

1 

1 

2  | 

1 

14.52891256 

1 

1 

3s2 . 3p5 . (2P*<3/2>) . 5p 

1 

1 

1 

2 [1/2] 

1 

1 

1 

1 

0  I 

1 

14.57594785 

1 

1 

1 

3s2 . 3p5 . (2P*<l/2>) . 5p 

1 

1 

1 

2 [3/2] 

1 

1 

1 

l 

1  1 

1 

14.68064940 

1 

1 

1 

3s2 . 3p5 . (2P*<l/2>) . 5p 

1 

1 

1 

2 [1/2] 

1 

1 

1 

1 

1  1 

1 

14.68711731 

1 

1 

1 

3s2 . 3p5 . (2P*<l/2>) . 5p 

1 

1 

1 

2  [3/2] 

1 

1 

1 

1 

2  | 
l 

14.68828937 

1 

1 

1 

3s2 . 3p5 . (2P*<3/2>) .4d 

1 

1 

2  [1/2]* 

1 

1 

1 

0  I 

14.6936388 

1 

1 

1 

1 

1 

1 

1  1 

1 

14.71089717 

1 

1 

3s2 . 3p5 . (2P*<l/2>) . 5p 

1 

1 

1 

2 [1/2] 

1 

1 

1 

1 

0  1 

1 

14.73811442 

1 

1 

1 

3s2.3p5. (2P*<3/2>) .4d 

1 

1 

1 

2  [3/2]* 

1 

1 

1 

1 

2  | 

1 

14.74253992 

1 

1 

1 

3s2.3p5. (2P*<3/2>) .4d 

1 

1 

2  [7/2]* 

1 

1 

1 

4  | 

14.7570507 

1 

1 

1 

1 

1 

1 

3  | 

1 

14.780511 

1 

1 

3s2.3p5. (2P*<3/2>) .4d 

1 

1 

2  [5/2]* 

1 

1 

1 

2  | 

14.809271 

1 

1 

1 

1 

1 

1 

3  | 

1 

14.824297 

1 

1 

3s2 . 3p5 . (2P*<3/2>) . 6s 

1 

1 

2  [3/2]* 

1 

1 

1 

2  | 

14.83881006 

1 

1 

1 

1 

1 

1 

1  1 
l 

14.84836805 

1 

1 

3s2 . 3p5 . (2P*<3/2>) .4d 

1 

1 

1 

2  [3/2]* 

1 

1 

1 

1 

1  1 
l 

14.859229 

1 

1 

1 

3s2 . 3p5 . (2P*<3/2>) .4f 

1 

1 

2 [3/2] 

1 

1 

1 

1  1 

14.9014407 

1 

1 

1 

1 

2  | 

14.9014901 

1 

3s2 . 3p5 . (2P*<3/2>) .4f 

1 

2 [9/2] 

5  | 

14.9037842 

1 

1 

4  | 

1 

14.9037893 

3s2 . 3p5 . (2P*<3/2>) . 4f 

1 

1 

2 [5/2] 

1 

3  I 

14.9065882 

1 

1 

2  | 
l 

14.9066270 

3s2 . 3p5 . (2P*<3/2>) . 4f 

1 

1 

2 [7/2] 

1 

3  I 

14.9090835 

1 

1 

4  | 

1 

14.9090883 

3s2 . 3p5 . (2P*<l/2>) . 4d 

1 

1 

1 

2  [3/2]* 

1 

2  | 

1 

14.95260379 

3s2 . 3p5 . (2P*<l/2>) . 4d 

1 

1 

2  [5/2]* 

1 

2  | 

14.95485109 

1 

1 

3  | 

l 

14.97152141 

3s2 . 3p5 . (2P*<l/2>) . 4d 

1 

1 

1 

2  [3/2]* 

1 

1  1 

1 

15.00356492 

3s2 . 3p5 . (2P*<3/2>) . 6p 

1 

1 

1 

2 [1/2] 

1 

1  1 

1 

15.0106092 

3s2 . 3p5 . (2P*<l/2>) . 6s 

1 

1 

2  [1/2]* 

1 

0  I 

15.0140645 

1 

1 

1  | 

1 

15.02208738 

3s2 . 3p5 . (2P*<3/2>) . 6p 

1 

1 

2 [5/2] 

1 

3  I 

15.0225900 

1 

1 

2  I 

1 

15.0258678 

3s2 . 3p5 . (2P*<3/2>) . 6p 

1 

1 

2 [3/2] 

1 

1  | 

15.0339710 

1 

1 

2  | 
l 

15.0356393 

3s2 . 3p5 . (2P*<3/2>) . 6p 

1 

1 

1 

2 [1/2] 

1 

0  I 

1 

15.0603891 

3s2 . 3p5 . (2P*<l/2>) . 4f 

1 

1 

2 [7/2] 

1 

3  I 

15.0830782 

1 

1 

4  | 

l 

15.0830834 

3s2 . 3p5 . (2P*<l/2>) . 4f 

1 

1 

2 [5/2] 

1 

3  I 

15.0832017 

1 

1 

2  I 

1 

15.0832446 

3s2 . 3p5 . (2P*<3/2>) . 5d 

1 

1 

2  [1/2]* 

1 

0  I 

15.10054275 

1 

1 

1  | 

1 

15.11774542 

3s2 . 3p5 . (2P*<3/2>) . 5d 

1 

1 

1 

2  [7/2]* 

1 

4  | 

l 

15.13054341 

3s2 . 3p5 . (2P*<3/2>) . 5d 

1 

1 

1 

2  [3/2]* 

1 

2  | 

1 

15.13684751 

3s2 . 3p5 . (2P*<3/2>) . 5d 

1 

1 

1 

2  [7/2]* 

1 

3  I 

1 

15.14592734 

3s2 . 3p5 . (2P*<3/2>) . 5d 

1 

1 

2  [5/2]* 

1 

2  | 

15.1610477 

1 

1 

3  I 

l 

15.1669454 

3s2 . 3p5 . (2P*<3/2>) .7s 

1 

1 

2  [3/2]* 

1 

2  | 

15.18063262 

1 

1 

1  1 

1 

15.18550704 

3s2 . 3p5 . (2P*<3/2>) . 5d 

1 

1 

1 

2  [3/2]* 

1 

1  1 

1 

15.1898229 

3s2 . 3p5 . (2P*<l/2>) . 6p 

1 

1 

1 

2 [1/2] 

1 

1  1 
l 

15.2006136 

3s2 . 3p5 . (2P*<l/2>) . 6p 

1 

1 

2 [3/2] 

1 

1  | 

15.2016566 

1 

1 

2  | 
l 

15.2048079 

3s2 . 3p5 . (2P*<3/2>) . 5f 

1 

1 

2 [3/2] 

1 

1  | 

15.2111325 

1 

1 

2  | 

1 

15.2111752 

3s2 . 3p5 . (2P*<3/2>) . 5f 

1 

1 

2 [9/2] 

1 

5  | 

15.2123191 

1 

1 

4  | 

1 

15.2123267 

3s2 . 3p5 . (2P*<3/2>) . 5f 

1 

1 

2 [5/2] 

1 

3  I 

15.2138399 

1 

1 

2  | 
l 

15.2138677 

3s2 . 3p5 . (2P*<3/2>) . 5g 

1 

1 

2  [5/2]* 

1 

3  I 

15.214032 

1 

1 

2  | 
l 

15.214032 

3s2 . 3p5 . (2P*<3/2>) . 5g 

1 

1 

2 [11/2] * 

1 

5  | 

15.214452 

1 

1 

6  I 

1 

15.214452 

3s2 . 3p5 . (2P*<3/2>) . 5f 

1 

1 

2 [7/2] 

1 

3  I 

15.2150700 

1 

1 

4  I 

i 

15.2150737 

3s2 . 3p5 . (2P*<3/2>) . 5g 

1 

1 

2  [7/2]* 

1 

4  | 

15.215274 

1 

1 

3  I 

1 

15.215274 

3s2 . 3p5 . (2P*<3/2>) . 5g 

1 

1 

2  [9/2]* 

1 

5  | 

15.215711 

1 

1 

4  | 

1 

15.215711 

3s2.3p5. (2P*<l/2>) .6p 

1 

1 

1 

2 [1/2] 

1 

0  I 

l 

15.2240866 

3s2 . 3p5 . (2P*<3/2>) .7p 

1 

1 

1 

2  [1/2] 

1 

1  1 

1 

15.2713708 

3s2 . 3p5 . (2P*<3/2>) .7p 

1 

1 

2 [5/2] 

1 

3  I 

15.2755689 

1 

1 

2  I 

l 

15.277418 

3s2 . 3p5 . (2P*<3/2>) . 7p 

1 

1 

2 [3/2] 

1 

1  | 

15.281660 

1 

1 

2  | 

1 

15.2824804 

3s2 . 3p5 . (2P*<l/2>) . 5d 

1 

1 

2  [3/2]* 

1 

2  | 

15.29629228 

Listing  1:  Partial  listing  of  ASCII  file  for  NIST  lines  of  Ar  I.  The  files  for  the  levels  have 
been  obtained  from  the  NIST  Atomic  Spectra  Database  (ASD)  website. 
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7  DATABASE  IMPLEMENTATION 


7.1.1  Electron-specific  levels 

The  initial  task  of  Colorado  when  first  executed,  is  to  read  in  the  electronic  levels  from 
an  ASCII  files  which  provides  the  configuration  name,  term,  J  value,  and  energy.  The 
partial  contents  of  such  a  file  is  provided  in  listing  7.1.1.  I/O  support  in  this  case  is 
provided  by  the  LevelReader  class  that  reads  in  the  file,  and  for  each  level,  instantiates  a 
new  Level  object.  The  reader  returns  a  Map  which  contains  a  collection  of  Level  objects, 
mapped  to  their  configuration  names,  facilitating  access  to  each  individual  level.  The 
class  Level  parses  and  stores  the  pertinent  information  including  configuration  names, 
quantum  numbers,  and  energies.  A  method  that  takes  a  collection  of  individual  levels 
and  combines  them  into  a  single  averaged  state  is  also  provided.  Partial  code  for  the 
LevelReader  and  Level  classes  are  provided  in  listings  1  and  2,  respectively. 

The  ASCII  files  containing  the  NIST  data  for  the  levels  are  named  according  to  the 
element  name  followed  by  Roman  numerals  for  the  ionization  level  with  the  file  name 
extension  .  levels.  For  example,  the  files  Arl .  levels,  Aril .  levels,  and  Arlll .  levels 
correspond  to  ground  state,  singly-,  and  doubly-ionized  argon,  respectively. 

Listing  1:  LevelReader  class 

public  class  LevelReader 

{ 

public  static  Map<S  t  r  ing  ,  Level>  read(String  specie  ){  ...  } 

} 


Listing  2:  Level  class 

public  class  Level 

{ 

/**  Unicode  descriptor  */ 
public  String  name; 

/**  NIST— assigned  configuration  */ 
public  S  tring  config  ; 

/ **  NIST— assigned  term  */ 
public  S  tring  term  ; 

/**  Valence  shell  *  / 
public  String  n£ ; 

/  **  Index  of  this  level  *  / 
public  int  index; 

/**  Principal  quantum  number  */ 
public  int  n  ; 

/**  Orbital  angular  momentum  quantum  number  */ 
public  int  £  ; 

/**  Total  orbital  angular  momentum  quantum  number  */ 
public  int  L; 

/*  *  */ 

public  int  K; 

/**  Total  angular  momentum  quantum  number  */ 
public  double  J  ; 

/  **  Statistical  weight  */ 
public  double  g  ; 

/  **  Statistical  weight  of  ionized  level  *  / 
public  double  gp  ; 

/ **  Number  of  optical  electrons  */ 
public  double  £ ; 

/ **  Excitation  energy  in  eV  */ 
public  double  e  ; 

/**  Ionization  potential  */ 
public  double  IP  ; 

/**  Core  quantum  number  */ 
public  double  jc; 

/**  Primed  (jc=l/2)  or  not  (jc=3/2)  */ 
public  boolean  primed  ; 

public  boolean  combined  ; 
public  Level  []  group; 

public  Level  (){} 

public  static  Level  parseNI  ST  (  S  t  r  i  n  g  config  ,  String  term,  double  J,  double  e){  ...  } 

public  static  String  getName  (  String  config  ,  String  term  ,  int  J){  ...  } 

public  static  Level  combine  (String  name  ,  Map<S  t  ring,  Level>  levels  ){  ...  } 

public  static  String  t  r  a  n  s  i  t  i  o  n  _t  y  p  e  (  Le  ve  1  £1  ,  Level  £2){  ...  } 

public  static  Level  getCopy  (  Level  level  ){  ...  } 
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7  DATABASE  IMPLEMENTATION 


7.1.2  Line  transitions 

In  much  the  same  way  as  a  collection  of  Level  objects  is  created  from  an  ASCII  file  of 
NIST  levels,  a  collection  of  Line  objects  is  created  from  a  file  of  NIST  transitions  lines. 
This  is  accomplished  by  the  LineReader  class  given  in  listing  3.  It  takes  as  an  argument 
the  collection  of  Level  objects  so  that  when  a  Line  object  is  instantiated,  it  includes  the 
necessary  information  about  the  transition  such  as  the  wavelength  as  well  as  pointers  to 
the  actual  upper  and  lower  Level  objects.  This  facilitates  access  to  the  level  information 
when  performing  operations  involving  the  transitions.  As  sample  fisting  of  the  NIST 
transition  lines  is  provided  in  listing  7.1.2  and  the  form  of  the  Line  class  is  given  in  listing 
4. 


Listing  3:  LineReader  class 

public  class  LineReader 
{ 

public  static  Map<  String  ,  Line>  read  (  String  specie  ,  MapK  String,  Level>  levels  ){  ...  } 

} 


Listing  4:  Line  class 

public  class  Line 
{ 

/**  Transition  wavelength  [m]  */ 

public  double  A; 

/  **  Einstein  coefficient  */ 
public  double  A; 

/**  Oscillator  strength  */ 
public  double  f  ; 

/**  Lower  level  */ 
public  Level  lower  ; 

/**  Upper  level  */ 
public  Level  upper  ; 

public  Line(){} 

public  static  Map<S  t  r  ing  ,  Line>  parseNIST(  String  path,  Map<S  t  r  i  ng  ,  Level>  levels  ){  ...  } 

public  static  Line  getCopy(Line  line){  ...  } 

} 
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7  DATABASE  IMPLEMENTATION 


00  00  00  00  00  CO  H  11)  Ifl  CO  N  H  H  11)  CO  ID  H  CD  ID  P)  CD  H  CO  ID  CO  ID  CO  CO  ID  CO  CO  ID  CO  H  S  CO  rl  ID  CO  ID  ID  CO  CO  CO  CO  H 

t — I  ■' — 1  i — 1  ■' — 1  ■' — 1  t— 1  CD  CD  ID  rl  ID  CO  00  CD  CO  CO  CO  CD  CD  CD  rl  CD  CO  ID  ID  CO  CD  CD  ID  ID  H  ID  rl  CD  ID  ID  CD  CD  CD  CD  CD  CO  CO  CD  rl  CD 


rH  rH  rH  rH  rH  -rH  O  Oi  M  rl  CD 

I  I  I  I  I  I  I  I  I  I  I 

OOOOO  O  H  rl  CN  O  N 


O  O  CN  rH  CN  O  1 


O  rl  Ol  H  N 
I  I  I  I  I 

H  rl  CN  CN  H 


CN  rl  O  CD  i 


I  CN  CN  O  CN  O  H  CN  CN 


O  CN  TH  CN  CN 
I  I  I  I  I 


C.O.O.C.C,  CN  CNCNCNCNCN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN  CN 

CO  H  CO  CO  rH  CO  H  CO  CO  CO  LO  H  H  CO  CO  LO  rH  rH  CO  rH  rH  rH  rH  CO  H  CO  rH  CO  CO  CO  rH  LO  CO  rH  LO  rH  H  CO  CO  LO  CO  rH  CO  rH  CO  rH 

CNCNCNCNCN  CN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN 


I  I  I  I 


I  I  I  I  I 


I  I  I  I  I 


3  3  3  3  S3  3  CN  CN  CN  CN  CN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN 

32222  22222  22222  ^  5  <2  S  S 

CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN  CNCNCNCNCN 


u  * 


8  5  m  &  $  $  t^COCOCOCO  S^IOIOID  COCOCOCOCO  CO  CO  LO  !0  !0  COCOCOLOLO  WCOCOCOLO  LOlOiOiOCO  COIOLOLOCO 

SSSSS  A  SSSSS  aaaaa  aaaaa  SSSSS  aaaaa  aaaaa  SSSSS  aaaaa 


m  %  nm  asm  nm  nm  m%% 

»»»»»  »  sssss  %%%%%  333”  *****  %m%  *****  sssss 


i  i  i  i  i 


i  i  i  I 


i  i  i  i  i 


Hill  |  55555  55555  55555  55555  55555  55555  55555  55555 

I I I I I  I  $$$$$  $$$$$  ^aaaa  aaaaa  aaaaa  aaaaa  aaaaa  aaaaa 

mu  uni  mu  uni  nm  nm  um  mu 

£££££  ggggg 

CO  co  co  co  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  co  co  co  co  cococococo  cococococo  cococococo  cococococo 

33333  3  3  3  3  3  33333  33333  55355  33333  33  3  3  3  3  3  3  3  3 


I  I  I  I 


Issli 


8  8  8  8  8  § 

OOOOO  O  j  j  j  j  ' 


:s  ss 

II II 

!  21  2  21  21  2!  2!  2! 


!  2  21  21  2  2!  2!  : 


12  2: 


OOOOO  O  QQQQQ  QQQQQ  QQQQQ  QQOOO  OOOOO  OOOOO  OOOOO  OOQQO 


iiiii  i  mu  fuii  iiiu  mu  nii!  mu  nm  mu 


i  CN  CO  ^  CN  CN  < 


LO  rH  rH  CO  rH  H  ID  (D  ■ 


I  CO  rH  CO  rl  d  H  CN  rl  00  (D  (O  H 


a-: 


Ill!!  I  mil  mil  Sill!  mil  Hill  IIIII  IIIII  Hill 

CD  CD  CN  N  ID  rH  CD  ID  CN  d  H  ID  S  CO  H  IN  OO^fCO^C~  N  ID  ID  N1  CO  CN  H  CN  ID  CN  ID  CN  O)  H  CD  CO  SCD  ID  CD  CN  H  N  H 

is-il  I  gsggg  nil  1  IIIII  IIIII  IIIII  IIIII 


ill 


IIIII  I  IIIII  IIIII  IIIII  IIIII  IIIII  IIIII  IIIII  IIIII 

»  — I  1  11333  33333  33333  IIIII  11555  55555  55555  55555 
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Listing  2:  Partial  listing  of  ASCII  file  for  NIST  lines  of  Ar  I. 


7  DATABASE  IMPLEMENTATION 


7.2  Input 


Once  the  collections  of  levels  and  lines  have  been  formed,  an  input  file  provided  by  the 
user  specifies  the  levels  of  particular  interest.  A  sample  input  file  for  Arl  is  provided  in 
listing  7.2.  This  particular  file  tells  ‘Colorado  to  use  only  the  first  31  levels  to  construct 
the  database.  Groups  of  levels  can  also  be  defined  by  including  incomplete  configurations 
such  as  4p  or  5  which  would  find  and  combine  all  available  levels  with  principle  quantum 
numbers  n  —  4  and  t  —  1  for  the  former  and  n  —  5  for  the  latter.  The  reader  together 
with  the  files  are  associated  with  the  namespace  Colorado  .  database .  nist .  levels.  The 
ordering  of  the  levels  in  the  input  file  is  important  since  it  determines  the  ordering  of  the 
rates  output  by  the  code.  It  is  not  a  requirement  for  the  levels  to  be  energy  ordered  or 
term  ordered.  The  input  file  is  Unicode  in  that  it  allows  sub-  and  superscripts. 


3p6 

4s[3/2]2 
4s[3/2]i 
4s’[1/2]0 
4s’  [1/2]  1 
4p[l/2]i 
4p[5/2]3 
4p[5/2]2 
4p[3/2] ! 

4p[3/2]2 
4p[l/2]0 
4p’[3/2]i 
4p’[3/2]2 
4p’[l/2]i 
4P’[l/2]o 
3d[l/2]0 
3d[l/2] ! 

3d[3/2]2 
3d[7/2]4 
3d[7/2]3 
3d[5/2]2 
5s[3/2]2 
5s[3/2]i 
3d[5/2]3 
3d[3/2] ! 

3d’[5/2]2 
3d’[3/2]2 
3d’[5/2]3 
5s’[1/2]0 
3s’  [1/2]  1 
3d’ [3/2]  i 

Listing  3:  A  sample  input  file  for  levels  of  Ar  I.  The  file  allows  the  user  to  communicate 
to  Colorado  which  levels  to  include  in  the  CR  database.  Unicode  characters  allow  clear 
configuration  descriptions  and  are  easily  handled  by  Java. 


7.3  Rate  calculations 

Calculation  of  reaction  rates  from  the  cross  sections  has  been  implemented  in  Colorado 
based  on  Maxwellian  distributions  of  both  the  heavy  particle  and  electron  energies.  With 
the  list  of  user-defined  levels,  the  code  computes  the  necessary  rates  for  each  of  the  possible 
processes  and  then  stores  them  as  a  function  of  temperature  in  separate  files,  numbered 
according  to  their  position  in  the  user-defined  input  file. 

7.4  Orgnaization 

The  cross  section  data  is  currently  (and  quite  tentatively)  organized  into  the  following 
namespaces.  Files  such  as  Arl.  levels  and  Arl.  lines  are  maintained  within  the  col- 
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orado.  database,  nist.  levels  and  Colorado  .database  .nist .  lines  namespaces.  User 
input  files  are  located  in  the  Colorado,  conf  ig  namespace  such  as  the  example  included 
in  listing  7.2.  The  actual  cross  sections  are  located  according  to  the  element  name 
and  process  type.  For  specific  models,  the  namespaces  are  labeled  accordingly  such 
as  Colorado . database . argon . deutsch  and  Colorado . database . argon . zatsarinny  for 
electron-impact  ionization  model  of  Deutsch  et  al  and  electron-impact  excitation  model 
of  Zatsarinny  &  Bartschat.  While  these  models  are  the  most  accurate  included  in  the 
database,  transitions  for  only  lower-lying  levels  are  available  and  therefore  Drawin’s  cross 
sections  are  used  to  complete  the  missing  information.  This  is  handled  in  Colorado  by 
first  looking  into  the  locations  of  the  specific  models  to  see  if  data  is  available,  and  then 
switching  to  the  general  models  when  the  data  cannot  be  found. 


Colorado 
Colorado . conf ig 

Colorado . database . argon . atomic_excitation 

Colorado . database . argon . atomic_ionization 

Colorado . database . argon . deutsch 

Colorado . database . argon . elastic 

Colorado . database . argon . electronic_excitat ion 

Colorado . database . argon . electronic_ionizat ion 

Colorado . database . argon . radiativerecombination 

Colorado . database . argon . zatsarinny 

Colorado . database . argon . zatsarinny_osc 

Colorado . database . krypton . atomic_excitation 

Colorado . database . krypton . atomic_ionization 

Colorado . database . krypton . deutsch 

Colorado . database . krypton . elastic 

Colorado . database . krypton . electronic_excitation 

Colorado . database . krypton . zatsarinny 

Colorado . database . nist . levels 

Colorado . database . nist . lines 

Colorado . database . xenon . atomic_excitation 

Colorado . database . xenon . atomic_ionization 

Colorado . database . xenon . deutsch 

Colorado . database . xenon . elastic 


Listing  4:  Organizational  structure  of  the  code  in  terms  of  namespaces. 


8  Conclusions  and  future  work 

Collsional- Radiative  models  for  the  noble  gases  have  been  detailed  which  included  atom- 
and  electron-impact  processes  as  well  as  radiative  processes.  These  models  are  intended  to 
support  plasmadynamic  simulations  under  conditions  of  high  non-equilibrium.  The  data 
has  been  presented  in  terms  of  angle-integrated  cross  sections  where  we  have  no  assump¬ 
tion  of  the  distribution  function  has  been  assumed.  This  leaves  open  the  possibility  to 
compute  the  CR  rates  using  an  arbitrary  distribution  function,  including  a  bin  approach. 

While  detailed  validation  has  been  provided  for  argon,  results  for  krypton  and  xenon  are 
still  preliminary.  Work  is  currently  underway  to  increase  the  level  of  sophistication.  These 
models  will  be  refined  and  benchmarked  against  further  experimental  data,  including 
expansion  flows  in  which  radiative  processes  dominate. 
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